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Abstract 

We consider the dynamics of a simple one dimensional model and 
we discuss the phenomenon of aging (i.e., the strong dependence of the 
dynamical correlation functions over the waiting time). Our model is 
the so-called random random walk, the toy model of a directed polymer 
evolving in a random medium. 
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Aging is a very interesting phenomenon that has been observed in spin 
glasses [1, 2, 3, 4], but it likely to be present in many other materials (for 
example in rubber), and to be a crucial signature of the behavior of a strongly 
disordered system. In brief we can say we have aging if the response of the 
system to a perturbation strongly depends on the time t w (waiting time) 
during which the system has been kept in the low temperature phase before 
starting the measurement. 

More precisely in an aging experiment one brings the system in the low 
temperature phase at time t = 0. One leaves the system quiescent up to 
the time t w , and at t w she adds a small perturbation (small enough that 
the linear response theory can be applied). One eventually measures the 
response function R(t, t w ) at time t + t w . Alternatively one can measure the 
correlation function C(t,t w ) among two quantities measured at time t w and 
at time t + t w . 

In the region where the observation time t is very smaller than the waiting 
time t w , t <^ t w , R and C are independent from t w and they are related by 
the fluctuation-dissipation theorem. Aging is present if R and C are strongly 
dependent on t w in the region where t is of the order of t w . The phenomenon 
of aging has been observed non only in real spin glasses, but also in numerical 
simulations of short range [5] and long range [6] spin glasses. In this last case 
the mean field approach [7, 8] is exact, so that an analytic study of aging 
should be possible. 

It seems very natural that in the region of large times a simple scaling 
behavior could hold. Different suggestions have been put forward. One 
simple possibility [9] is that in the region where both t and t w are large 

C(t,t w )=t-»f{±-) , (1) 

where the exponent // may be equal to zero. Unfortunately at the moment 
we do not have a complete justification of this formula. 

Waving our hands we can sketch a simple qualitative explanation. We 
suppose that the system is described by a potential with a corrugated land- 
scape, where many barriers of all possible heights are present. Such a system 
(considered in the thermodynamic limit) is always out of equilibrium. Its 
dynamics is given by the never ending search of the absolute minimum by 
crossing potential barriers that can be arbitrary large. 

Let us assume that the available phase space does not increase too much 
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with the height of the barriers that have been crossed. In this case we can 
suppose that at time t w and temperature T the system has explored all the 
phase space that can be reached from the origin by crossing barriers that 
are smaller than T\n(t w ). The system therefore will remain near the bottom 
of the explored phase space up to the moment at which it crosses one other 
large barrier. This can happen only at times of the order of t w . If we assume 
that the shape of the deepest minimum found does not depend on the size 
of the explored region we are lead to eq. (1). In this way we can explain 
the dependence on j-. One would need a more careful analysis to derive 
the value of fi, which may depend on the detailed quantity whose correlation 
function is computed. It should be noted that if 

lim lim C(t, t w ) = const ^ (2) 

t — >oo tw — 

(i.e. if C goes to a non zero finite limit in the region where both t and t w are 
large and t <C t w ) then ji = 0. 

The aim of this note is to discuss the behavior of a very simple system, 
the one dimensional random walk in random environment [11]. In this simple 
case many analytic results are available (for a review see [10]). Our numerical 
simulations will show that the aging phenomenon is present also in such a 
simple setting. 

The model we consider is a particle which equilibrium probability distri- 
bution is proportional to 

P( X ) OC e-^M+A* 2 ) ; (3) 

where V is a random Gaussian quantity with zero average and correlations 

(V{x)-V{y)f = \x-y\ a . (4) 

In the case a = 1 (the one that we will mainly consider in the following) 
the force F(x) = — ^ is uncorrelated from point to point. One introduce 
the contribution \x 2 to regularize the static equilibrium distribution, and 
A — > in the relevant limit. At equilibrium this model coincide with the 
toy model that was introduced in the study of the behavior of interfaces in 
random medium (see [12] and references therein). 

It is known that [13] 
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(x 2 ) - (x) 2 ~ A" 1 , 
ln((:r 2 ) - (x) 2 ) ~ 0(1). 



(5) 



The first equation implies that for non-zero A the particle is contained 
inside a region of radius of order A~ 2 / 3 . In the A — > limit the particle is 
always at infinity at equilibrium (the potential V is unbound). 

In most cases the particle is localized around the minimum. The most 
probable value of {x 2 ) — (x) 2 is of order one, although rare events, in which the 
potential has two widely separated nearly degenerate minima, do contribute 
strongly to the average of (x 2 ) — (x) 2 , making it of order A -1 . 

The dynamics of this model is a random walk biased to produce the 
equilibrium distribution (3) at large times. More precisely the model can be 
defined by introducing a one dimensional lattice. If the particle is at point 
x at time t (x and t taking integer values) then p x is the probability that at 
time t+1 the particle is at site x + 1, and 1 — p x the probability that it is 
at site x — 1. p x is a random variable with zero average (typically uniformly 
distributed in (0, 1), not including the interval limits). 

It is easy to see that the equilibrium potential corresponding to this dy- 
namics satisfy the condition 

V{x + 1) - V(x) = - ln( Px ) . (6) 

1 — Px+l 

Consequently at large distance the difference of the potential V(x) — V(y) 
are approximately Gaussian, with a — 1. 

This dynamics has been widely investigated [10]. It can be proven [14] 
that the particle to go from x to y has to cross (with high probability) a 
barrier that is of order \x — y\ 1 ^ 2 . If the system is in x = at t = 0, at a 
large time t it will be in x ~ ln(t) 2 . 

In the following we will discuss the phenomenon of aging in this model. 
The simplest quantity one could study is the correlation function 

C(t, t w ) = (x(t + t w )-x(t w )) 2 , (7) 

that for large t (in the asymptotic region studied by Sinai, where t ^> t w ) 
behaves as ln(t) 4 . But this is not a good choice. In the region of large times, 
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for t <C t w , this function behaves as ln(t) 3 . The reasons for this behavior are 
quite clear. If we assume that the probability distribution for the dynamics 
at time t may be mimed by the probability distribution for the statics, where 
A is chosen so that the value of (x) coincides with the correct one, we find 
from the previous equations that (a; 2 ) — (x) 2 is of order ln(t) 3 . In other words 
the system has explored a region of size ln(t) 2 and the probability of having 
two nearly equal minima inside this region vanishes follows from 

the analysis of ref. [13]. 

The previous analysis implies that the aging properties of the function C 
are not simple and the scaling law in eq. (1) cannot hold with /i — 0. So we 
found convenient to consider the following correlation 

L(t, t w ) = \n((x(t + t w ) - x(t w )) 2 ) , (8) 

or equivalently M(t,t w ) = exp(L(£, t w )). M is not far from the most likely 
value of (x(t + t w ) — x(t w )) 2 . In this case we expect that in the region of large 
times, when t <^t w , 

M(t, t w ) = constant + 0(m(£) _1 ) . (9) 

All that said, it seemed natural to set up a numerical simulation to verify 
the validity of the scaling law: 

M(t,t w ) = f(^) + 0(\n(t)- 1 ) . (10) 

We have generated 15 x 10 3 realizations of the one-dimensional random 
potential. For each of these realizations of the random potential we have 
observed the random walker travelling for 2 K steps and we have measured 
the correlation M(t, t w ) at times t and t w equal to 2 fc , with k=\,...,K—\. 
The spatial lattice had an infinite extent, i.e. the particle would never hit a 
boundary. We have taken K = 21, i.e., performed order of 2 x 10 6 sweeps 
per sample. 

We show in Fig. 1 (ln(|x(i)|)) versus ln(ln(t + 1)) for different waiting 
time t w , scaled logarithmically. In Fig. 1 we have drawn continuous lines 
joining the points of equal t w . Higher curves correspond to lower t w . The 
first continuous line from the top joins points taken after no waiting (t w = 0), 
the second the ones with t w — 1, the third the ones with t w — 2, the fourth 
the ones with t w — 4 and so on, with an exponential progression. 
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Figure 1: The average logarithm of the distance reached by the walker after 
t steps as a function of ln(ln(t)). Here the lines join points with constant t w . 
Lines from top to bottom are from lower to higher t w . 
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The upper lines of Fig. 1 have t t w , and in the right part of the plot 
they scale with good precision as expected. Such straight lines have indeed 
a slope very close to 2. 

In Fig. 2 we plot the same points, but we join points of constant j- ratio. 
Here the j- scaling is quite clear. Indeed the lines at intermediate values 
of j- tend for large t to a constant value. The distance reached from the 

tw 

walker at time t only depends over the ratio j-. If we double the observation 

tw 

time t but we also double the waiting time t w the average distance covered 
by the walker does not change. Figures 1 and 2, and the other figures we 
will show, allow us to draw the main conclusion of this note. For large t and 
t w , t of order t w , the correlation functions behave as an universal function of 
j-. The simplicity of the Id (toy) model allowed us to get a numerical very 
compelling evidence for such an effect. 

In Fig. 3 we plot again the same data, but this time we select the ones with 
both t and t w larger than 32, and we plot exp((ln(|x(t)|))) versus ln(^±^-). 
The data fall with good precision on a single scaling curve. In the left part 
of the (small t as compared to t w ) the universal curve is a constant 

(aside from small corrections), while in the large j- region it behaves as 
ln(^) 2 (according to the Sinai [14] scaling law). 

Similar results can be obtained considering the energy-energy correlation 
defined by 

C E (t,t w ) = (V(x(t)) -V(x(t w ))f , (11) 

where V is defined by equation (6). 

We expect that CE(t,t w ) ~ ln(^) 2 , independently on the value of a. For 
the same reason we also expect that the expectation value of V: 

E(t)=WWfi) (12) 

behaves as — ln(t) for large t, independently on a. 

Such a scaling law is very well satisfied. In Figures 4 — 6 we give the 
analogous of Figures 1 — 3, and we exhibit the analogous scaling laws (which 
in this case are even better than in the former one). In Fig. 4 we give —V(t) 
as a function of ln(i), and we join points of equal t w . The large t linear 
behavior is very clear. In Fig. 5 the same points, but we join point of equal 
j- ratio. Again, the curves tend to constant values. In Fig. 6 the rescaling of 

tw 
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Figure 2: Same points of Fig. 1, but here the lines join points with constant 
_t_ 
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log((t+l)/(t w +l)) 

Figure 3: A selection of the points of Fig. 1 and 2, for large t and t W) versus 
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the points with t > 32, t w > 32, which again work with very good accuracy. 
Here it is very clear, for example, that in the t < t w asymptotic region the 
walker does not gain any energy. 

Let us conclude. The numerical results we have presented here bring good 
evidence for the correctness of the simple aging scaling relation 

C(t,t w )~F(^) , (13) 

for the correlation function of the logarithm of the distance and of the energy. 

In our model simple aging is correct. This conclusion seems to us interest- 
ing because of its potential general implications, and as far as the simplicity 
of the model is such that a more sound analytic (and may be rigorous) deriva- 
tion of these results does not seem impossible. 

As a further development of this work, we notice it could be interesting to 
study finite volume effects for large, but finite lattice size (i.e., with reflecting 
or periodic boundary conditions) and to study the modification of our results 
for very large time. The extension of the model to higher dimensional cases, 
where the structure of minima is more complex, will be probably instructive. 



10 



15 - 



10 - 

A 
■+-> 



5 - 



- 




5 10 

log(t+l) 

Figure 4: As in Fig. 1 (lines join points with constant t w ), but for —V(t) 
versus \n(t). 
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Figure 6: As in Fig. 3, but for — V(t). 
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